Revisiting the Hanbury Brown-Twiss Setup for Phase Fluctuating Bose Gases 
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The Hanbury Brown-Twiss experiment has proved to be an effective means of measuring two- 
point correlation function of identical particles. We analyze experimental observation of stripes 
formation of a phase fluctuating Bose- Einstein condensates in a highly elongated 3D traps [Dettmer 
et al,, Phys. Rev. Lett. 87, 160406 (2001)] by means of axial two-point correlation functions. We 
also predict that the stripes are present in quasi-lD Bose gas in the mean-field as well as in the 
hard-core bosons regimes. These stripes can be realized by measuring the axial two-point correlation 
functions by using the Bragg interferometric method which is similar to the original Hanbury Brown 
and Twiss experimental setup. 

PACS numbers: 03.75.Hh, 03.75.Nt 



I. INTRODUCTION 



Since the pioneering works on the realization of Bose- 
Einstein condensates (BEC) [1] of alkali- atoms, a great 
variety of experimental and theoretical investigations 
have probed the macroscopic phase coherence of confined 
quantum gases. For a trapped 3D BEC well below the 
transition temperature T c , experiments have confirmed 
the macroscopic phase coherence by measuring the cor- 
relation length which is equal to the condensate size [2]. 
However, the phase coherence strongly depends on the 
shape of the confining potential which can be control at 
will. It was theoretically proposed that the axial phase 
fluctuations of an elongated 3D BEC can be very large 
in the equilibrium state, where the density fluctuations 
are strongly suppressed [3]. The axial phase coherence 
length in the elongated systems can be smaller than the 
axial size and this is referred as the quasi-condensates. 
Similarly, quasi-lD Bose gas in the mean-field regime [4] 
forms a quasi-condensate [5] with a large phase fluctua- 
tions even at temperatures as low as 0.1T C [6]. The quasi- 
lD Bose gas behaves like a hard-core bosons (commonly 
known as Tonks- Girardeau Gas) [7] when the strength 
of the two-body potential is very strong which has been 
observed experimentally [8]. The hard-core bosonic sys- 
tems have a very strong vacuum phase fluctuations which 
prevents from being a (quasi-)condensate even at T = 
and mimicking the exclusion principle for fermions. The 
most striking features of the phase fluctuating 3D BEC 
is that the stripes formation in the axial density profile 
of the system [9] . 

The celebrated Hanbury Brown and Twiss (HBT) ex- 
periment [10] in which the spatial second-order correla- 
tion function 62(2) of a light source is characterized by 
measuring the correlations of intensity fluctuations in the 
wave field. The original idea of the HBT experiment is 
that by measuring the intensity correlations between two 
separated beams, they essentially compared the inten- 
sities at two different points in the unseperated beam 



[11]. The equal-times second-order correlation function 
provides information on the relative spatial distribution 
of pairs of identical particles. Therefore, the stripes for- 
mation can be understand by analyzing the second-order 
correlation function if it shows oscillatory behavior. 

In this Letter, we understand the experimental obser- 
vations of the stripes formation in the axial density profile 
of a phase fluctuating 3D BEC by means of the axial two- 
particle correlation functions. Next, we propose that the 
same stripes formation in the axial density profile can be 
observed in quasi- ID Bose gas in the mean- field as well as 
in the hard-core bosons regimes. We also discuss how to 
realize the stripes formation by a suitable choices of the 
Bragg pulse in the Bragg interferometer which is similar 
to the HBT experiment. 

Two-point correlation function: We consider a 3D 
BEC confined in an elongated harmonic trap, where 
the repulsive mean-field interaction energy exceeds the 
radial (TkJo) an d the axial (hio z ) trap energies. At 
T = 0, the density profile has the well-known form 
n (p,z) = (fi/g)(l - P 2 /Rl - z 2 /Zl), where (j, = 
0.5huj z [15Nau^/a z ujl] 2 ^ 5 is the zero-temperature chem- 
ical potential, and g = Anah /m is the two-body in- 
teraction strength. Also, Rq = y / 2/i / 'muj^ and Zq = 
•y/2/x/mLjJ are the radial and the axial size of the con- 
densate, respectively. Due to the repulsive mean-field 
energy, density fluctuations are strongly suppressed in 
a trapped BEC. Therefore, the bosonic field operator 
describing the condensate can be written in the form 
ijj(r) = \/no (r)exp[i(f>(r)], where the phase operator <f)(r) 
is defined as 
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Here, d„ is the annihilation operator of the quasiparti- 
cle excitation with quantum numbers v and energy hio v \ 
ip v is the corresponding quasiparticles normalized wave 
functions. 

The normalized two-particle correlation function is de- 



fined as 

4 

C2[{r,:}] = niVM^r 1 < ^(ri)^(r 2 )V>(r 3 Mr4) > 

2 = 1 

Due to the strong suppression of the density fluctua- 
tions, the normalized density correlation function of the 
trapped condensate is constant, i.e. C2(ri,r2,r2,ri) = 1. 
Therefore, a simple measurement of the density corre- 
lations in the condensate is not enough to describe co- 
herence properties. However, by measuring the density 
correlations in the interference pattern generated by two 
spatially displaced copies of a parent BEC, it is possi- 
ble to correlate the bosonic field operator «/>(r) at four 
different positions and extract the C^ri,^,^,^) [12]. 

3D cigar-shaped BEC: As described in Ref. [9] , the 3D 
condensate consists of TV = 5 x 10 5 atoms of 87 Rb. The 
radial and axial trapping frequencies are ujq = 2-7TX 365 Hz 
and lo z — 27rx 14 Hz (a z — ^Jh/mLo z = 2.879/im), respec- 
tively. The condensate is already elongated along the lon- 
gitudinal z axis due to the aspect ratio A = u>q/uj z = 26. 
In addition, further radial compression or expansion of 
the condensate are obtained by applying the superim- 
posed blue detuned optical dipole trap to the magnetic 
trap. In the experiment [9], they performed the mea- 
surements for the radial trapping frequencies ujq between 
2ir x 138 Hz and 2ir x 715 Hz corresponding to aspect 
ratios A between 9.8 and 51, respectively. 

The excitations of a cigar-shaped BEC can be divided 
into two regimes: "low-energy" axial excitations with en- 
ergy Hlu z < E a << huJo, and "high-energy" radial excita- 
tions with E r > hu>o >> huj z . As pointed out in Ref. [3] 
that the low-energy axial excitations have wave-lengths 
larger than Rq and exhibit a pronounced ID behavior 
which gives the most important contribution to the low- 
energy axial phase fluctuations. The low-energy axial 
modes have the energy spectrum e/ = fouj z y/j(j + 3)/4 
[13]. The normalized wave functions ipj of these quasi- 
particle modes have the form 
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where P- ' (z) are Jacobi polynomials and z = z/Zo is 
a dimensionless variable. 

Using Eq. (1) for the phase operator, F({ri}) becomes, 
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hujj 



U ZoRl (j + l)V7(7T3)"^ L 2fc S T J 
x [PJ^fr) + PJ^fr) 

-(lf' 1 \zs)+Pi 1 ' 1 \S4))] a . (4) 

The upper cut-off limit (j maK ) on the summation can 
be obtained from the constraint: ej < [i. It is use- 
ful to define the space variables in the following way: 



z.\ = (s + d)/2 — — Z2, and zj, = (s — d)/2 = —S4. Here, 
d is the displacement between the two interfering con- 
densate copies and s is the relative distance between the 
positions in the interference pattern at which the particle 
densities are evaluated. This new choice of variables has 
been realized in the experiment [12]. 

In Fig.l and Fig. 2, we plot the axial two-particle cor- 
relation function vs the relative separation at T = and 
T = 0.6T C , respectively, for three different choices of the 
aspect ratio. 
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FIG. 1. (Color online) Plots of the axial two-particle cor- 
relation functions vs the relative separation. 
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FIG. 2. (Color online) Plots of the axial two-particle cor- 
relation functions vs the relative separation. 



The axial two-particle correlation function at zero tem- 
perature has oscillatory behavior with very small ampli- 
tude. It shows little evidence of the stripes formation 
even at zero temperature. Nevertheless, we do not claim 
that stripes are present even at T = because it is almost 
one for a large separation. The minima of the C^s) at 
various relative separations are decreasing with the in- 
creasing of the aspect ratio A. It implies that vacuum 
phase fluctuation increases with the aspect ratio. 




Fig. 2 shows the oscillatory behavior with large am- 
plitudes compared to the zero temperature case. These 
prominent deep valleys implies the presence of stripes in 
the axial density distribution due to the strong thermally 
excited phase fluctuations and the large aspect ratio A. 
The minima of the C*2(d, s) at various relative separations 
are decreasing with the increasing of the temperature and 
the aspect ratio A. 

In the actual experiment [9], the stripes formation is 
observed after the ballistic expansion of the atomic cloud. 
However, it does not rule out the possibility of the stripes 
formation in static BEC i.e., without ballistic expansion. 
According to the Ref . [9] , the observation of the stripes 
formation after the ballistic expansion of the cloud is 
due to the rapidly decreasing of the mean-field interac- 
tion and the axial velocity fields arc then converted into 
the density distribution. In our studies we find that the 
stripes formation in the axial density profile is already 
present even before switching off the trap. Therefore, 
the stripes formation are not necessarily due to the axial 
velocity fields which are converted into density modula- 
tions during the ballistic expansion, as stated in Ref. [9]. 
The width of the stripes are too small to be observed 
directly in-trap due to the lack of the experimental reso- 
lution, but it is seen after expansion of the cloud since the 
stripes in the static cloud is enlarged during the ballistic 
expansion of the cloud. 

Quasi-ID Bose gas in the mean-field regime: Now we 
consider quasi- ID Bose gas in the mean-field regime as 
described in Ref. [4]. In particular, the system consists 
of N ~ 10 4 atoms of 23 Na in the trap with axial trap- 
ping frequency u) z = 2n x 3.5 Hz, and radial trapping 
frequency luq = 2n x 360 Hz. The low-energy exci- 
tation spectrum is given by luj — ui z \/j{j + l)/2 and 
the corresponding normalized cigcnfunctions are 4>j{z) — 

J^±fPj(z/Z ), [14] where Pj(z/Z ) is the Legendre 

polynomial in z and Zq — a z (2.998Naa z /aQ) 1 ^ 3 is the 
Thomas-Fermi half length. 

Using Eq. (1) for the phase operator, F({ri}) for 
quasi-lD system becomes, 

x \Ptfx) + P,{h) - Pj(~z 3 ) - Pjih)} 2 . (5) 

Using the expression for F({zj}), we plot two-point cor- 
relation functions vs the relative distance at T = and 
T = 0.6T C in Fig. 3 and Fig. 4, respectively. 

Similar to the 3D BEC in very elongated traps, the 
vacuum phase fluctuations are not able to produce the 
stripes in quasi-lD Bose gas in the mean-field regime, 
However, at large temperature, the two-point correla- 
tion functions have prominent oscillatory behavior with 
large amplitudes, compared to the zero temperature case. 



Therefore, the stripes are present in the quasi-lD Bose 
gas in the mean-field regime at large temperature. 
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FIG. 3. (Color online) Plots of two-point correlation func- 
tions vs the relative distance for various trapping aspect ratio. 
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FIG. 4. (Color online) Plots of two-point correlation func- 
tions vs the relative distance for various trapping aspect ratio. 

Quasi-lD Bose gas in the hard-core bosons regime: Let 
us consider quasi-lD Bose gas in the hard-core bosons 
regime. The low-energy excitation spectrum is given 
by ej = jhuj z and the corresponding normalized eigen- 

functions are ipj(z) = ,/ -^-Tj(z/Zq), where Tj(z/Zo) is 
the first-order Chebyshev polynomial in z. Here, Zq = 
y/2Na z is the Thomas- Fermi half-length of the hard-core 
bosons system with the chemical potential /i = Nhuj z 
[15]. 

Using Eq. (1) for the phase operator, F"({ri}) for quasi- 
lD system in the hard-core bosons regime at T = be- 
comes 

N 1 
F[ Zl ] = J2 -P^i) + Tj{h) - T 3 (z 3 ) - T 3 (z 4 )} 2 . (6) 

3=1 J 

Fig. 5 shows the two-point correlation function of the 
quasi- ID Bose gas in the hard-core regime at T = 0. 
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FIG. 5. (Color online) Plots of two-point correlation func- 
tions vs the relative distance for various number of hard-core 
atoms. 

Note that Fig. 5 is valid only when s > 1/2N due to 
the hydrodynamic approximation. This hydrodynamic 
approximation is failed to describe the short-range cor- 
relations due to the strong interactions when s < 1/2N. 
Fig. 5 shows that the two-point correlation functions 
have the oscillatory behavior with large amplitude even 
at zero temperature. This implies that stripes are al- 
ready present in the hard-core bosons regime even at zero 
temperature due to the strong quantum phase fluctua- 
tions. Note that the amplitude of the oscillations in the 
zero temperature two-point correlation functions is very 
very small in the cigar-shaped 3D BEC as well as in the 
quasi-lD Bose gas in the mean field regime, compared 
to the quasi-lD Bose gas in the hard-core regime. The 
two-point correlation functions with large amplitude of 
the oscillations implies the fermionic like behavior in the 
hard-core Bose gas at zero temperature. 

Detection: The realization of the stripes formation is 
possible by measuring the axial two-particle correlation 
function. The Bragg interferometric method presented in 
Ref. [12] to measure the two-particle correlation function 
is analogous to the original HBT experiment [10]. There- 
fore, the prominent valleys in the axial two-particle cor- 
relation function should be observe by using the Bragg 
interferometric method as described in Ref. [12] with a 
proper choices of the wave vector k and the time inter- 
val (At) between the two Bragg pulses such that the 
distances (d = 2TikAt/m) between two auto correlated 
copies should be d ~ 0.1 — 0.2. We have checked that 
for other choices of d, d > 0.2, the oscillatory behavior 
in the two-point correlation functions washes out. This 
is due to the fact that the two interfering matter waves 
are not HBT-correlated when d > 0.2. Here we mean 
HBT-correlated in a sense that the matter waves will 
produce an HBT effect at the output ports, as opposed 
to the different question of the correlations in the parent 
condensate [11]. In fact, the two matter waves are HBT- 



correlated only if they are interfering within a coherence 
time (r), the characteristic time scale in the experiment. 
The relative distance d between two interfering BECs at 
the output ports of the interferometer can be controlled 
by varying the time interval (At). When At is small 
such that d ~ 0.1 — 0.2, we expect the correlations to be 
maximal and shows the HBT effect. For large At > r, 
the distance between two interfering BECs is large so 
that it becomes HBT-uncorrelated and then it does not 
produce valleys in the interference patterns. The time in- 
terval, At, before colliding two copies should be less than 
the coherence time r which is inversely proportional to 
the temperature of the condensate i.e. At < t ~ h/hsT, 
where T is the temperature of the parent condensate [11]. 
In conclusion, we have understand the stripes in the ax- 
ial density profile of a strong phase fluctuating 3D cigar- 
shaped BECs in terms of the two-point correlation func- 
tion. We have also predicted that the stripes are present 
in the quasi-lD Bose gas in the mean-field regime at large 
temperature as well as in the ID hard-core bosonic sys- 
tems even at zero temperature due to to the strong quan- 
tum phase fluctuations. We have also pointed out that 
one can realize the stripes by a proper choices of the pa- 
rameters in the Bragg interferometer experiment which 
is similar to the HBT setup. This experiment would pro- 
vide the direct realization of the stripes formation, which 
is of fundamental importance in the phase fluctuating 
BECs. 
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